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ABSTRACT The fungal pathogen Cochliobolus carbonum (anamorph, Bipolahs zeicola) causes Northern 
Leaf Spot, leading to a ubiquitous and devastating foliar disease of corn in Yunnan Province, China. Asexual 
spores (conidia) play a major role in both epidemics and pathogenesis of Northern Leaf Spot, but the 
molecular mechanism of conidiation in C. carbonum has remained elusive. Here, using a map-based cloning 
strategy, we cloned a single dominant gene, designated as BZconl (for Bipolahs zeicola conidiation), which 
encodes a predicted unknown protein containing 402 amino acids, with two common conserved SANT/ 
Myb domains in N-terminal. The BZconl knockout mutant completely lost the capability to produce con- 
idiophores and conidia but displayed no effect on hyphal growth and sexual reproduction. The introduced 
BZconl gene fully complemented the BZconl null mutation, restoring the capability for sporulation. These 
data suggested that the BZconl gene is essential for the conidiation of C. carbonum. 
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In many phytopathogenic fungi, asexual sporulation in the form of 
conidia is critical in the life cycle, and so deciphering the molecular 
mechanisms involved in conidiation development is a prerequisite to 
provide novel strategies for disease management. The survival of these 
fungi is ensured by the production and dissemination of spores that 
are long-lived and highly resistant to environmental stresses. Spore 
germination represents a critical stage in the life cycle of fungi and 
constitutes a prerequisite for colonization in a new environment 
(Lafon et al. 2005). The process of conidiation involves many com- 
mon developmental themes, including gene expression, cell speciali- 
zation, and intercellular communication. For example, in Aspergillus 
nidulans conidiation, three central regulatory genes (brlA, abaA, and 
wet A) control conidiation- specific gene expression and determine the 
order of gene activation during conidiophore development and co- 
nidial maturation (Park and Yu 2012). These studies suggest that it is 
important to develop new strategies to control the occurrence of 
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fungal disease by deciphering the molecular mechanism involved in 
asexual sporulation of plant pathogenic fungi. 

Asexual spores play a major role in the epidemics of Cochliobolus 
carbonum (anamorph, Bipolaris zeicola), which causes Northern Leaf 
Spot, a ubiquitous and devastating foliar disease of corn in many 
regions of the world (Nelson 1959, 1970; Welz and Leonard 1994; 
Zhang et al. 2013). In the wild, C. carbonum exists in the asexual form, 
reproducing mainly through the production of conidia. As soon as the 
conidia are in contact with the corn surface, they release a copious, 
two -layered, extracellular matrix (ECM) resulting in the adhesion of 
germlings of the fungus to the host surface (Evans et al. 1982; Apoga 
et al. 2001). After penetration of the germ tubes, infectious hyphae 
initially grow within the plant cells. Eventually, the plant cells are 
killed and thousands of conidia could be released from each lesion. 
Despite the importance of conidia in pathogenic epidemics, little is 
known about the molecular mechanism of conidiation in C. carbonum. 

In ascomycetes, conidia formation requires the temporal and 
spatial control of cell differentiation, which is a process under 
polygenic control (Adams et al. 1998; Etxebeste et al. 2010). To date, 
no genes involved in conidiation have been reported in the asexual 
development of C. carbonum. Therefore, it would be valuable to iden- 
tify and characterize genes that are specifically involved in conidia 
formation for an in-depth understanding of the asexual development. 
Development of genetic maps and the identification of molecular 
markers closely linked to target genes have allowed the map-based 
cloning of specificity genes governing the phenotype (Sweigard et al. 
1995; Farman and Leong 1998; Ma et al. 2006). Here, via both a 
map -based cloning strategy and functional experiments, this study 
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characterizes a novel C. carbonum gene, BZconl, which is a determi- 
nant of conidiation in this fungus. BZconl encodes a potential binding 
protein containing 402 amino acids as putative transcription factors 
and regulates the conidiation in C. carbonum. 

MATERIALS AND METHODS 

Fungal strains, cultures, and transformation 

Field isolates S92 and SI 29 were collected from corn plants with 
Northern Leaf Spot symptoms in Yunnan Province, China, in 2010. 
Experimental strains Aa82, Aal08, and Aall3 were obtained from the 
progeny of the cross S92 X SI 29. Other strains in this study are shown 
in Table 1. Mating types of all isolates were determined by polymerase 
chain reaction (PCR) amplification as described previously by Gafur 
et al. (1997). Mycelia collected from potato dextrose agar (PDA) plates 
were used for isolating fungal genomic DNA and RNA. Mycelia col- 
lected from Fries medium (Fries 1978) were used for obtaining the 
protoplasts. Protoplasts were isolated and transformed as described by 
Liu and Friesen (2012). The selection medium was supplemented with 
100 |xg/ml hygromycin B (Roche, USA) or 500 |mg/ml G418 (Amer- 
esco, USA) for the screening for hygromycin B-resistant (HR) or 
G4 18 -resistant (GR) transformants, respectively. Single conidium of 
the transformant was isolated to ensure nuclear homogeneity (Choi 
et al. 1999). If lack of conidia precluded single-spore isolation, then 
transformants were transferred three to four times to select against 
heterokaryons (Lev et al. 1999). Isolates and mutants were stored on 
filter paper disks as previously described (Silue et al. 1992). 

Cross and conidiation assay 

The isolates S92, Aa82, and Aall3 of C. carbonum produced abun- 
dant conidia (con Y ), whereas the isolate SI 29 did not produce conidia 
(con N ) on PDA plates. Crossing was performed according to the de- 
scription of Welz and Leonard (1994) on Sachs agar medium (Gafur 
et al. 1998). Pseudothecia were harvested 20 to 25 days after crossing, 
and single ascospores were isolated from the asci and allowed to 
germinate on PDA. For conidiation assays, all isolates were incubated 
on PDA plates (6-cm diameter) at 25° for 10 d in the dark. Then, the 
colony surface was washed with 10 ml sterile distilled water (SDW) 
and the conidia suspension was filtered through two layers of cheese- 
cloth. The amount of conidia per milliliter of the suspension was 
determined by counting the conidia with a hemacytometer using light 
microscopy. An isolate was regarded as a con Y or a con N progeny with 
the presence or absence of conidia, respectively. Each isolate was 
tested at least four times. Segregation of con Y in the progeny popula- 
tion was examined with a x 2 test for goodness of fit for con Y /con N 
ratios. A summary of the segregation populations is shown in Table 2. 



Linkage marker screening of the con Y locus 

Pools of DNA for bulked segregant analysis (BSA) (Michelmore et al. 
1991) were prepared by combining equal amounts of DNA from five 
con Y or five con N progeny strains. RAPD (Welsh and McClelland 
1990; Williams et al. 1990), ISSR (Zietkiewicz et al. 1994), and SRAP 
(Li and Quiros 2001) amplifications were performed on each pool and 
parental strains. RAPD reaction productions were analyzed by elec- 
trophoresis on agarose gels (1.2%). The specific fragment for con Y 
pools was excised from the agarose gels. ISSR and SRAP reaction 
productions were resolved in polyacrylamide gels (6%) coupled with 
silver- staining. The exclusive band for con N pools was excised from 
polyacrylamide gels as described by Zhang et al. (2009). The PCR 
production was excised and purified using the MinElute Gel Extrac- 
tion Kit (Sangon Biotech, China). The eluted fragment was ligated into 
pMD18-T simple vector (Takara, Japan) following the supplier's 
instructions and transformed into competent Escherichia coli TGI, 
and then sequenced by Beijing Genomics Institute (BGI; Beijing, 
China). The sequence obtained was used to design primers for se- 
quence characterized amplified regions (SCARs) (Paran and Michelmore 
1993) analysis. The successful SCAR primer sequences used are shown 
in Table 3. 

Genetic mapping of the gene BZconl 

For genetic linkage analysis, SCAR amplifications were performed 
with genomic DNA from 271 progenies of the cross between Aall3 
and SI 29. Segregation data obtained were analyzed with the software 
Mapmaker version 3 (Lander et al. 1987). Parameters for map con- 
struction were a minimum logarithm of odds (LOD) score of 3.0 and 
a maximum recombination fraction of 0.4. The Kosambi function was 
used to compute the recombination distances in centimorgans (cM). 

Construction and screening of BAC library 

The Aall3 genomic DNA bacterial artificial chromosome (BAC) 
library was constructed as described by Brosch et al. (1998) 
with minor modifications. After Hindlll partial digestion, DNA frag- 
ments 75 to 120 kb in size were ligated into plasmid pBeloBACll and 
electroporated into E. coli TGI using an electroporator (Eppendorf, 
Germany). After electroporation, the cells were resuspended in 600 jixl 
of LB medium (10 g tryptone, 5 g yeast extract, and 10 g NaCl in 1000 
ml deionized water), allowed to recover for 1 h at 37° with gentle 
shaking, and then plated on LB agar containing chloramphenicol (12.5 
|jLg/ml). The plates were incubated overnight, and the recombinants 
were picked manually and transferred to 96- well plates that were in- 
cubated for 10 h with shaking at 220 rpm. The BAC library was stored 
at —80°. Screen of the BAC clone by pool method was first performed 



Table 1 Wild-type and mutant strains of Cochliobolus carbonum used in this study 



Strains 



Brief Description 



Growth (cm) a 



Conidiation (x10 4 ) b 



S92 

S129 

Aa82 

Aa108 

Aa113 

M113-1 

C113-2 

C129-1 



Field isolate, MAT1-2-1 

Field isolate, MAT1-1-1 

Experimental isolate, MAT1-2-1 

Experimental isolate, MAT1-1-1 

Experimental isolate, MAT1-2-1 

BZconl knockout mutant of Aa1 13 

BZconl complemented transformant of M 1 1 3-1 

BZconl complemented transformant of S129 



7.52 : 

7.55 

7.62 

7.62 

7.53 

7.56 

7.51 

7.54 



0.10 c 
: 0.21 
: 0.04 
: 0.03 
: 0.12 
: 0.02 
: 0.06 
: 0.02 



101.34 ± 0.0 C 



0 



120.12 
110.11 
110.23 



0 



12.13 
10.25 



0.04 
0.06 
0.03 

0.04 
0.06 



Vegetative growth was measured at day 6 after incubation on PDA plate. 

Conidiation was measured by counting the number of conidia collected with 1 0 ml of sterilized distilled water from 6-cm diam. PDA plate at day 1 0 after incubation 
at 25°. 

Data were presented as mean and SDs were calculated from four independent experiments. 
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Table 2 Segregation of random ascospore progeny from crosses between conidia-producing isolates (con Y ) and isolates producing no 
conidia (con N ) of Cochliobolus carbonum 



Cross Code 


Cross (con Y xcon N ) 


Generation 


No. of Progeny 
Con Y Con N 


Expected Ratio 


X 2 


Probability 


Aa 


S92xS129 


F1 


158 


150 


1:1 


0.20 


0.65 


BC 


Aa82xS129 


Backcross 


26 


24 


1:1 


0.08 


0.78 


DD 


Aa113xS129 


Backcross 


134 


137 


1:1 


0.03 


0.86 


KC 


Aa108xM113-1 




59 


55 


1:1 


0.14 


0.71 



with the markers SC403 and SC254. Gaps present in the assembled 
sequence were further filled as described by Zhou et ol. (2006). The 
selected BAC clones were sequenced, and the obtained sequences were 
assembled using the Bioedit software (http://www.mbio.ncsu.edu/ 
bioedit/bioedit.html) . 

Chromosome walk and gene prediction 

We adopted chromosome walking of physical analyses to delimit the 
location of the BZconl gene. For the first move, we constructed BAC 
contig mapping and confined gene BZconl to a genomic interval 
between markers SC403 and SC254, and then this interval was ana- 
lyzed using SSR finder (http://www.fresnostate.edu/ssrfinder). For the 
second move, obtained SSR sites were used to design primers detect- 
ing polymorphism between the parents, and the screened primers 
were further used to amplify the special fragment of genomic DNA 
from 271 progenies of the cross between Aal 13 and SI 29, respectively. 
A marker, SSR2 (Table 3), developed based on SSR-PCR amplifica- 
tion, further delimited BZconl to a 54,526-bp genomic interval. Genes 
of the derived interval sequence were predicted using the Fgenesh 
program (http:/ /linuxl .softberry.com/berry.phtml) . 

Nucleic acid manipulation 

Total DNA of the isolates and mutants was extracted and purified 
from the mycelium by using the cetyltrimethylammonium bromide 
(CTAB) protocol (He 2000). PCR amplifications were performed in 
a total volume of 20 jjlI containing 0.4 |jlM of each dNTP, 5 fxM of 
each primer, 1 unit of easyTVzg or 2 units easyp/w DNA polymerase 
(Trans, China), 2.0 jjlI of 10X reaction buffer, and 10 to 20 ng of 
genomic DNA. Total RNA was extracted from the mycelia grown 
for 3 to 8 d on the PDA plates using the Trizol reagent (Invitrogen, 
USA) according to the manufacturer's protocol. Reverse-transcription 
PCR (RT-PCR) was performed with the PrimeScript strand cDNA 
Synthesis Kit (Takara, Japan) following the supplier's instructions. 
Procedures for gel electrophoresis, plasmid DNA preparations, and 
restriction enzyme digestions were performed following standard pro- 
tocols (Smith 1993; Stacey and Isaac 1994). For Southern blot analysis, 
genomic DNA was isolated from C. carbonum wild-type strain Aal 13 
and putative BZconl deletion mutant Ml 13-1. Aliquots of 15 mg 
DNA were digested with Sail, separated by electrophoresis in 0.8% 
agarose gels, and transferred onto a Hybond N + membrane (Amer- 
sham Pharmacia Biotech). The probe was amplified with primers 
H-F2 and H-R2 (Figure 2A and Table 3), whereas probe labeling and 
hybridization were performed with Amersham ECL Direct Nucleic Acid 
Labeling and Detection system (GE Healthcare UK Limited Companies), 
following the manufacturer's instructions. All primers (Table 3) were 
synthesized by BGI. 

Knockout of gene BZconl 

For BZconl knockout, the replacement vector (pBSKCl) favoring 
double-crossover integration was constructed in three steps (Figure 2A). 



First, a 1084-bp Xhol-Hindlll upstream fragment of the BZconl 
gene was ligated to X/zoI-HmdIII-digested pBlueScript II KS (+) 
(pBS), and the ligated plasmid was denominated as pBSPll. Then, 
a 1317-bp Hindlll-Pstl downstream fragment of the BZconl gene 
was ligated to the HmdIII-P5£l-digested pBSPll, and the ligated 
plasmid was designated as pBSPHP12. Subsequently, a 2.0-kb hygB 
resistance gene cassette (hph) y digested with Hindlll from plasmid 
pTFCM (Li et al 2005), was ligated to Hindlll- digested pBSPHP12. 
Homologous fragments used in the construction of the knockout 
vector were amplified using Aal 13 as the template with two pairs 
of primers, C-F1/C-R1 and C-F2/C-R2, respectively (Table 3). Cir- 
cular pBSKCl was transformed into Aal 13, and targeted integration 
of the plasmid into the C. carbonum genome was confirmed by DNA 
gel blot analysis. The transformants were further determined by PCR 
with the primer pairs CF1/CF2 (Table 3). 



Table 3 List of oligonucleotides used in this study 


1^ cim c 




SC254 


F, TG G GTCCCCAG CAAACCG 




R, TGGGTCCCCCTG CTT AG C 


SC383 


F, TGGGTCCCTCTTACC 1 1 1 1 G 




R, TG G G TCCCTC AC AATATG CC 


SC403 


F, GGGGGATGAGGCGGGTAGAG 




R, G G G G G ATG AG ACCG ATAACAGT 


SC829 


F, GTGTGTGTGTGTCGGATGG 




R, GTGTGTGTGTGGCCTGAATCAA 


SC855 


F, ACACACACACACCTG G CAC 




R, ACACACACACACACACCTCCTT 


SC24 


F, GAATTTGCCCGCA 1 1 1 1 I CC 




R, GTCCAAACCGGACCCTGAAC 


S8-35 


F, AAGCTTGGGTTCTTACGCGA 




R, CCTCGAATGGACCCGAC 


S8-33 


F, AAG CTTATCG G CG AACTG C 




R, CGTGATGGAGCACGCGA 


S1-125 


F, AAG CTTCTCCCG G AAG CTATC A 




R, G T AAG TTCTCCCTTC ACTC AAG ATT 


S1-123 


F, AAGCTTTGTCAATAGAATCC 1 1 1 1 1 GG 




R, CATCCCCTCACTTGTAAATACACAT 


SSR2 


F, TG CT ATTG TCGTCGCCGC 




R, ATTTTATTTTATTTTA 


CF1 


ATG GTTG CCCACAG ACG 


CR1 


CTACATGATGTTTCGCAGAG 


C-F1 


CCG CTCG AG AG AACG CTG 


C-R1 


AATAAG CTTG G G CCCTG AG CAG CAAAA 


C-F2 


AATAAG CTTAATAG CCG AGTTG TCCACCA 


C-R2 


AAAG AATTCTAG CAACG G CAT AG CAACA 


HB-CF1 


AAAGGATCCGCCGAAACGAGAACAAGATG 


HB-CR1 


AAAG G ATCCTCG AG AG AACG CTG CCT 


H-F2 


ATG AAAAAG CCTG AACTCAC 


H-R2 


CTATTCCTTTG CCCCTCG 


F. forward primer; 


R. reverse primer. 


The sequence u 


nderlined indicates the restriction cut site for construction of 


knockout and complementation vector. 
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Figure 1 Genetic and physical maps around the BZconl locus in the isolate Aa113. Genetic maps were constructed using multipoint linkage 
analysis on segregating populations of 271 progeny with the Mapmaker/Exp 3.0 software. Dashed lines connect the physical location of markers 
to their location on the genetic map or connect markers to their location in the bacterial artificial chromosome (BAC) clone contig. The gene was 
predicted according to Fgenesh programs and the arrows correspond to position of gene candidates. CAN13 is the gene identified as a candidate 
for BZconl and used for knockout and complementation assays. 



Complementation of gene BZconl 

For complementation assays, a BZconl region including the BZconl 
along with its 1.7-kb upstream and 0.4-kb downstream sequences was 
amplified using Aall3 as the template with primer pair HB-CF1/HB- 
CR1 (Table 3). It was subsequently cloned into the BamHl sites in 
plasmid pFA6a-kanMX6 harboring a G418 resistance cassette (KanMX) 
as pFACON3L Circular pFACON31 was transformed into the BZconl 
mutant M113-1. 

RESULTS 

Map-based cloning of the gene BZconl 

The progeny resulting from three independent crosses segregating for 
conidiation was summarized in Table 2. The segregation of conidial 
and nonconidial strains in the Fl progeny (cross S92/S129) and BC 
populations (Aa82/S129 and Aall3/S129) fitted a ratio of 1:1 ratio in 
the x 2 test of goodness-of-fit, indicating that the asexual sporulation of 
C. carbonum is controlled by a single dominant gene. We designated 
this gene as BZconl (for Bipolar is zeicola conidiation). 

A PCR-based marker system was used to detect polymorphism 
between two parental isolates. Using a BC family of 271 strains 
produced by crossing the con Y isolate Aall3 and con N isolate S129, 
a total of six SCAR primers were found to be linked to the targeted 
BZconl gene. By using software MapMaker 3.0, the BZconl locus was 
mapped between SC403 and SC254, which were identified as 0.6 and 
3.0 cM, respectively (Figure 1). 

The construction of a BAC contig including the BZconl locus in 
con Y parental isolate Aall3 was initiated by SC403 and SC254 



markers. By PCR screening, SC403 and SC254 marker fragments were 
detected in BAC clones BG8-3 and BA1-12, respectively. So, BG8-3 
and BA1-12 were completely sequenced and a gap was present be- 
tween BG8-3 and BA1-12. The BG8-3 and BA1-12 end sequences 
were then used to design several contiguous markers (Figure 1 and 
Table 3) and allowed the identification of a BAC clone (HD6-1) 
bridging the gap. HD6-1 was also completely sequenced. As a result, 
the gene BZconl interval, between markers SC403 and SC254, con- 
sisted of three overlapping BAC contigs, which covered a distance of 
217,367 bp (Figure 1), giving an average physical-to-genetic-distance 
ratio of 60.4 kb/cM. 

Forty- four SSR sites were obtained from the BZconl interval 
sequence using SSR finder. The SSR-PCR amplifications showed 
that only marker SSR2 was identified as co-segregation with 
BZconl. By using MapMaker 3.0 software, the BZconl locus was 
determined, located further between markers SC403 and SSR2, 
and the physical distance between SC403 and SSR2 was 54,526 
bp. A total of 13 genes, designated CAN1 to CAN 13, were pre- 
dicted in this region (Figure 1). 

These genes were used to design gene code-specific primers (data 
not shown) and applied to the amplification of con N isolate S129. The 
obtained sequences were determined. It was found that the coding 
sequences in the con Y isolate Aall3 were identical to the ones in con N 
isolate S129, indicating that the transcriptions of genes may be differ- 
ent. RT-PCR experiments showed CAN2, CAN10, and CAN12 were 
not detected in both of these two strains, and CAN1, CAN3, CAN4, 
CAN5, CAN6, CAN7, CANS, CAN% and CAN 11 were both present in 
Aall3 and SI 29. The remaining CAN 13 was conspicuously absent in 
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C-Fl 
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C-Rl 
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CR1 C-F2 
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C-R2 
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11 
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B 
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H-R2 
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7.1 kb 
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Figure 2 Strategy for BZconl gene knockout and con- 
firmation of BZconl mutant. (A) BZconl knockout vec- 
tor pBSKO contains a 2.0-kb hygB resistance gene 
cassette. Homologous recombination through a double 
crossover event results in the knockout of a BZconl re- 
gion with the hygB resistance cassette (hph). This re- 
gion includes a 634-bp-long upstream fragment, the 
entire open reading frame (1 209-bp-long) of BZconl 
gene, and a 76-bp-long downstream fragment. The 
positions of primers are indicated by small arrows. X, 
Xhol. H, HindW. E, EcoR\. (B) Southern blot analysis. 
Total genomic DNA samples isolated from Aa113 (ex- 
perimental isolate) and M 1 1 3-1 {BZconl knockout mu- 
tant of Aa1 13) were digested with Sa/I. The hph prober, 
a 1 -kb PCR fragment amplified from pBSKO using pri- 
mers H-F2 and H-R2, is exactly the BZconl region 
replaced by the hygB resistance gene. (C) Total RNA 
extracted from the mycelia of Aa113, M 11 3-1, and 
C1 13-2 grown for 3 to 8 d on PDA plates, respectively. 
And obtained RNA was subjected to RT-PCR using 
BZconl gene-specific primers CF1 and CR1 . The RT- 
PCR product is an approximately 1 .2-kb fragment in 
Aa113 and C113-1 {BZconl complemented transform- 
ant of M1 13-1) as predicted, but is missing in M1 13-1 . 



the con N isolate SI 29 in RT-PCR experiments, showing that CAN 13 
could be the only BZconl gene candidate. 

Determination of BZconl by gene knockout and 
complementation tests 

To assess whether CAN 13, as BZconl gene candidate, is responsible 
for the conidiation of the fungus, homologous recombination was 
used to replace CAN 13 with the plasmid pBSKCl, which contains 
the 2.0-kb hygB resistance cassette. The wild-type strain Aall3 was 
transformed with pBSKCl, and 26 HR transformants were isolated. 
All the transformants did not form conidia on the PDA. Genetic 
crosses were performed between HR transformants, the mutant 
Ml 13-1 and an experimental strain, Aal08 (Figure 3). We then iso- 
lated 1 14 single ascospore strains from the progeny of this cross. All 
55 HR strains failed to produce conidia, but the 59 hygromycin- 
sensitive (HS) strains were normal in conidia production. PCR am- 
plification with CF1 and CR1 revealed that all HS strains had an 
approximately 1.2-kb band that was absent in all HR strains. This 
proves to be only one locus replaced by the hygB resistance gene. 



Cosegregation of the replacement event and HygB resistance was also 
confirmed by Southern bolt analysis (Figure 2B). 

To further confirm that the disruption of the CAN 13 was respon- 
sible for the phenotype of the CAN 13 mutant, we conducted a com- 
plementation assay by transferring the vector pFACON31 harboring 
the full-length CAN 13 and a G4 18 -resistant gene as a selectable 
marker into the mutant, Ml 13-1. One G418-resistant complementary 
transformant, CI 13-2, fully complemented the CAN 13 null mutation, 
restoring the capability for sporulation (Table 1 and Figure 4B). Like- 
wise, the vector pFACON31 was introduced into the con N strain 
SI 29, and the transformant, C 129-1, gained the capability to produce 
conidia on PDA. Additionally, based on RT-PCR analysis, CAN 13 
transcripts were detected in Aall3 and CI 13-2, but were absent in 
Ml 13-1 (Figure 2C). These data clearly demonstrated that CAN13 and 
BZconl was the same gene. 

BZconl is required for conidiation 

After growing on PDA plates for 6 d, all colony diameters were 
determined by using at least four replicates for each strain or mutant. 



S92XS129 



Aall3XS129 



Aa82XS129 



Aal08XM113-l 




Figure 3 Mating result of the BZconl knocked-out mu- 
tant on Sachs agar medium. Pseudothecia is visible as 
dark spots. S92, field isolate, MAT1-2-1; S129, field iso- 
late, MAT7-7-7; Aa1 13, experimental isolate, MAT1-2-1; 
Aa82, experimental isolate, MAT1-2-1; Aa108, experi- 
mental isolate, MAT1-1-1; and M113-1, MAT1-2-1, 
BZcon 1 knockout mutant of Aa1 1 3. 
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A Aall3 M113-1 C113-2 S129 C129-1 




B Aall3 M113-1 C113-2 S129 C129-1 




Figure 4 Phenotypes of the wild-type strain Aa113 and S129, BZconl knockout mutant M 113-1, BZconl complemented transformants C113-2 
and C129-1. (A) Colony morphology of Aa113, S129, M113-1, C113-2, and C129-1 on potato dextrose agar (PDA) plates inoculated for 6 d. (B) 
Plate cultures observed under light microscope (red arrows showing conidia). 



The BZconl deletion mutant, Ml 13-1, displayed no obvious defects in 
the mycelial growth rate when compared with the strain Aal 13 (Table 
1), but Ml 13-1 produced gray- white aerial hyphae similar to that of 



the SI 29, in contrast to Aal 13, which had gray-green aerial hyphae 
(Figure 4A). In sporulation assays, Ml 13-1 completely lost the capa- 
bility to produce conidia on PDA. When the aerial hyphae were 



12? _bP_ ^ + 436 bp 

Aal 1 3 BZconl Allele T AAC T CC TCG ACCGCC TTGGCG TCC AAG TGTCAGGACTACTCCCTCCT CACGCCCGC CTTG AT ATATACCGGTGCTTCGCCTCT ACC ACGTTT ACCGTTGGTC CGTC AGACAACGTTC ACTACC ACTC AACATCCACT CG 

S129 BZconl Allele TCAACATCCACTCG 

BZconl CDS 

Aal 13 BZconl Allele ACTCAACAACACACTTGCACACACTTTGAATCTACCCTTCCACAGTCACCAACCCGCACATACGGTATGTCTCCTCTCTCTCCCTTCTTTAGACCCTTTGCACCATCTCGCACtGCCCCGGGTACATCCCGTCCTCTGTA 

S129 BZconl Allele ACTCAACAACACACTTGCACACACTTTGAATCTACCCTTCCACAGTCACCAACCCGCACATACGGTATGTCTCCTCTCTCTCCCTTCTTTAGACCCTTTGCACCATCTCGCACTGCCCCGGGTACATCCCGTCCTCTGTA 
BZconl CDS 

Aal 1 3 BZconl Allele TTTTGTCTGC AACAC ATCTTGTCAATCTAC ACATTGTC AGTGCCTGTTGC AGATAT ACCCCCTCCCTTCTCTC ATCTCCCCTGTTTCC ACCTATCAACGCTC AGCCTC ATGGTTGTCTCACCCCCCTCGCATTTTCATCT 

S129 BZconl Allele TTTTGTCTGC AACAC ATCTTGTCAATCTAC ACATTGTC AGTGCCTGTTGC AGATAT ACCCCCTCCCTTCTCTC ATCTCCCCTGTTTCC ACCTATCAACGCTC AGCCTC ATGGTTGTCTCACCCCCCTCGCATTTTCATCT 

BZconl CDS 

Aal 1 3 BZconl Allele TGTGTCATGGTTACCTCC AGTCCACGAAACTGACCTGCC ACCGCTTAGCCTTCCATCGAACATCGACTCAGAGTCAAAGTTCGCATTTCTTCTC ACTTGTTTTCAAAC AACTGC ACCCTTTCTTAACCACAAAACTTC AC 

812 ft BZconl Allele TGTGTCATGGTTACCTCCAGTCCACGAAACTGACCTGCCACCGCTTAGCCTTCCATCGAACATCGACTCAGAGTCAAAGTTCGCATTTCTTCTC ACTTGTTTTCAAAC AACTGC ACCCTTTCTTAACCACAAAACTTC AC 

BZconl CDS 

ATG— > 

Aal 1 3 BZconl Allele AA TG G TT GCCCAC AGA CGC GGC CC C TGG TC G CAG CACGAG GACC AGTGGCT TGT CAA CC TCGTAGC TC G CAATG GAC CC C ACAAC T GGG TTCGCATTT CACAAGAAATC GGCTCC AGG TC ACC CAAG CAA TG CAG AGA G C 

S 1 2 9 BZconl Allele AATGGTTGCCCAC AGACGCGGCCCCTGGTCGC AG CACGAG GACC AGTGGCTTGTCAACCTCGT AG C TCG CAATG GAC CCC ACAACTGGGTTCGC ATTTCACAAGAAATCGGCTCC AGGTC ACC CAAG CAATG CAG AGAGC 

BZconl CDS -ATGGTTGCCCACAGACGCGGCCCCTGGTCGCAGCACGAGGACCAGTGGCTTGTCAACCTCGTAGCTCGCAATGGACCCC ACAACTGGGTTCGC ATTTCACAAGAAATCGGCTCC AGG TCACCCAAGCAATGCAG AGAGC 



Aal 13 BZconl Allele 
S129 BZconl Allele 
BZconl CDS 



GCTATCACCAGAACCTCAAGCCCTCATTAAACCACGACCCAATCACTGGCGAAGAGGGCGAACTGATTGAGAAGATGGTTGCTGAGATGGGTAAGAGATGGGCAGAGATTGCTCGCCGGTTAAAAGGTCGCAGCGACAAT 
GCTATCACCAGAACCTCAAGCCCTCATTAAACCACGACCCAATCACTGGCGAAGAGGGCGAACTGATTGAGAAGATGGTTGCTGAGATGGGTAAGAGATGGGCAGAGATTGCTCGCCGGTTAAAAGGTCGCAGCGACAAT 
GCTATCACCAGAACCTCAAGCCCTCATTAAACCACGACCCAATCACTGGCGAAGAGGGCGAACTGATTGAGAAGATGGTTGCTGAGATGGGTAAGAGATGGGCAGAGATTGCTCGCCGGTTAAAAGGTCGCAGCGACAAT 



Aall3 BZconl Allele 
S129 BZconl Allele 
BZconl cos 



GCTGTCAAGAACTGGTGGAACGGAGGTCAAAACCGCAGAAAGCGCAATCCTGAGCGTGCCGACCAACGACCCCGCGACTACATGCACCAAGCACAGCCTGGCTTCCATTACCCACTACCTAGCACTCTTCCTCACCATGA 
GCTGTCAAGAACTGGTGGAACGGAGGTCAAAACCGCAGAAAGCGCAATCCTGAGCGTGCCGACCAACGACCCCGCGACTACATGCACCAAGCACAGCCTGGCTTCCATTACCCACTACCTAGCACTCTTCCTCACCATGA 
GCTGTCAAGAACTGGTGGAACGGAGGTCAAAACCGCAGAAAGCGCAATCCTGAGCGTGCCGACCAACGACCCCGCGACTACATGCACCAAGCACAGCCTGGCTTCCATTACCCACTACCTAGCACTCTTCCTCACCATGA 



Aal 13 BZconl Allele 
S129 BZconl Allele 
BZconl CDS 



AGAGCACCACCACATGACCCCTCGCCATGACCAGCAACGACACGACTCCCACCGACCGTCAGTACACTCGCCACATGCTCTCAAGCTTCCCCAGCCCCAGCGATTTGAGCCCTACCCAGCACCTCGCCCACAGCCTCTCG 
AGAGCACCACCACATGACCCCTCGCCATGACCAGCAACGACACGACTCCCACCGACCGTCAGTACACTCGCCACATGCTCTCAAGCTTCCCCAGCCCCAGCGATTTGAGCCCTACCCAGCACCTCGCCCACAGCCTCTCG 
AGAGCACCACCACATGACCCCTCGCCATGACCAGCAACGACACGACTCCCACCGACCGTCAGTACACTCGCCACATGCTCTCAAGCTTCCCCAGCCCCAGCGATTTGAGCCCTACCCAGCACCTCGCCCACAGCCTCTCG 



Aal 13 BZconl Allele 
S129 BZconl Allele 
BZconl CDS 



ATCTCGCCACATACGGATCACACGCGGCTGTACGAGCACGAGGATTCGAGACGCCTATGCCTTCACCTTCTGGATACTCCATCGCATCTGCCGACGGCGCTCCTTCCCTCATCACTGACACTGGCTCAGAGTCGAGGTCC 
ATCTCGCCACATACGGATCACACGCGGCTGTACGAGCACGAGGATTCGAGACGCCTATGCCTTCACCTTCTGGATACTCCATCGCATCTGCCGACGGCGCTCCTTCCCTCATCACTGACACTGGCTCAGAGTCGAGGTCC 
ATCTCGCCACATACGGATCACACGCGGCTGTACGAGCACGAGGATTCGAGACGCCTATGCCTTCACCTTCTGGATACTCCATCGCATCTGCCGACGGCGCTCCTTCCCTCATCACTGACACTGGCTCAGAGTCGAGGTCC 



Aall3 BZconl Allele 
S129 BZconl Allele 
BZeonl CDS 



CCCCGTGGTGCCGCTTCTCCCCCAGACATGACACTCGCACCGCTAGTCGGTAGCCGGGAAGAGCGCAAGAACTCCAGCACCCGTTACCTCCCCGCCACTGGCTTTGCACCTGAAGACGACTCCTCCAGCGATCCTTACCG 
CCCCGTGGTGCCGCTTCTCCCCCAGACATGACACTCGCACCGCTAGTCGGTAGCCGGGAAGAGCGCAAGAACTCCAGCACCCGTTACCTCCCCGCCACTGGCTTTGCACCTGAAGACGACTCCTCCAGCGATCCTTACCG 
CCCCGTGGTGCCGCTTCTCCCCCAGACATGACACTCGCACCGCTAGTCGGTAGCCGGGAAGAGCGCAAGAACTCCAGCACCCGTTACCTCCCCGCCACTGGCTTTGCACCTGAAGACGACTCCTCCAGCGATCCTTACCG 



Aal 13 BZconl Allele 
S129 BZconl Allele 
BZconl CDS 



ACGCCGTGGCTCACTTCAGCTGCCTCCTCTGCACCGCAACGAGGCCCAGACTTCCAAGGCAGCCCCTTCTTTCC GCGAGCCAAACTACGAG CGCCCGCAACCACCACAACTCCTTACCCAGCCGACCCCCCAGCTCACAC 
ACGCCGTGGCTCACTTCAGCTGCCTCCTCTGCACCGCAACGAGGCCCAGACTTCCAAGGCAGCCCCTTCTTTCCGCGAGCCAAACTACGAGCGCCCGCAACCACCACAACTCCTTACCCAGCCGACCCCCCAGCTCACAC 
AC GC C GT GGC T C ACTT C AG CTG CC T CC T CTG C AC CGC A AC G AGGC C C AG AC TTC CAA GGCAGCCCCTTCTTTCCGCGAGC C A AAC TACGAGCGCCCGCAACCACCACAACTCCTTACCCAGCCGACCCC CC AGCTC AC AC 



Aal 13 BZeonl Allele 
S129 BZconl Allele 
BZconl CDS 



ACCACCACTCTATGCCTCAGTCTGCTCCCACCAGTGTACCACAACGACCACTACCACCTTTTCACTCACTCGCTGAGCCAAGGCTACCCCTGCCCCACTCCACAAGCGCGAGCAGCAGCCCGGCATCGAGGCAACTGCCG 
ACCACCACTCTATGCCTCAGTCTGCTCCCACCAGTGTACCACAACGACCACTACCACCTTTTCACTCACTCGCTGAGCCAAGGCTACCCCTGCCCCACTCCACAAGCGCGAGCAGCAGCCCGGCATCGAGGCAACTGCCG 
ACCACCACTC T ATGCC TCAGTC TGCTCCCACCAGTGT ACC ACAAC GAC CAC TAC CAC CTT TT CAC T C AC TC GC T GAG CCA AGG CT A CC C CT G CC C CACTCC ACAAGCGC GAGCAGC AGCCCGGCATCGAGGCAACTGC CG 



Aal 13 BZeonl Allele 
S129 BZconl Allele 
BZeonl CDS 



TCCCCCTACGTCTTGCCCTCAATGCCAGCAGAGGCACGAAGCGCACCCGTATCGCCAAAGGACCGAATGTCTCTGCGAAACATCATGTAG 
TCCCCCTACGTCTTGCCCTCAATGCCAGCAGAGGCACGAAGCGCACCCGTATCGCCAAAGGACCGAATGTCTCTGCGAAACATCATGTAG 
TCCCCCTACGTCTTGCCCTCAATGCCAGCAGAGGCACGAAGCGCACCCGTATCGCCAAAGGACCGAATGTCTCTGCGAAACATCATGTAG 



Figure 5 Sequence comparisons of the upstream region of BZconl from two isolates A113 and S129. Dot line with double arrows shows the 
absence of 126 continuous nucleotides in the upstream (+426 bp) of BZconl from the isolate S129 when compared with that of Aa1 13. 
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Figure 6 An alignment of the amino acid sequences of BZCON1 and its 11 related homologs. Two SANT/Myb-like domains were located at 
N-terminal (approximately 100 amino acids in length). 



examined, abundant conidiophores and conidia were produced by 
Aall3, but none at all were produced by Ml 13-1 mutant. These 
results indicated that the BZconl was essential for conidiation of C. 
carbonum, and the development of conidia and conidiophores was 
blocked in the mutant (Figure 4B). 

A short upstream promoter region potentially mediates 
transcriptional regulation of the gene BZconl 

Approximately 1.7-kb upstream sequences of BZconl from the isolates 
SI 29 (not produce conidia in PDA; without transcription of BZconl) 
were amplified and sequenced. Sequence alignment found the absence 
of 126 continuous nucleotides in the upstream of BZconl from the 
isolate S129 when compared with that of Aall3 (Figure 5). We further 
observed that transformant C 129-1 can gain the ability to produce 
conidia on PDA via acquiring a segment including BZconl with the 



1.7-kb upstream sequence of Aall3. The results indicated that the 
deficiency of a short putative promoter 126 bp in length may contrib- 
ute to limitation of BZconl transcription. 

BZconl encodes a predicted unknown protein with two 
conserved SANT/Myb-like domains at N terminal 

The ORF of BZconl gene (GenBank accession number: KJ562713) 
was 1209 bp long, encoding 402 amino acids, and there was no intron 
in the ORF. To confirm the absence of introns, a full-length BZconl 
complementary DNA (cDNA) was obtained from Aall3 by RT-PCR 
using primers CF1/CR1. The amino acid sequence of BZconl was 
analyzed using InterProScan (Hunter et al. 2012), showing BZconl 
has two SANT/Myb domains in the N-terminal region (approximately 
100 amino acids in length). This result indicated that BZconl should 
be a putative Myb-like transcription factor. To further strengthen this 
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prediction, using amino acid sequence of BZCON1, we performed 
blast searches (E- value <1.00e-5) against nonredundant protein se- 
quence database (ftp://ftp.ncbi.nlm.nih.gov/blast/db/) and revealed 
that all 164 top -blast hits were derived from two related fungal phy- 
lum Ascomycota (162 hits) and Basidiomycota (2 hits) (Supporting 
Information, Table SI and Figure SI), having two common conserved 
SANT/Myb-like domains at N terminal (Figure S2). Surprisingly, the 
top nine blast-hits have similar amino acid sequence lengths (399-402 
aa) and 73%- 100% amino acid sequence identities as BZCON1 (Fig- 
ure 5 and Table SI). BZCON1 and its top nine blast-hits were clus- 
tered into a group with bootstrap support of 100%, belonging to three 
fungal family of subclass Pleosporomycetidae (Figure SI). However, 
the most related to the group lost approximately 200 to 300 amino 
acid residues in the C-terminal region (e.g., Tuber melanosporum and 
Pyronema omphalodes in Figure 6 and Figure SI). 

DISCUSSION 

Conidiogenesis is a widespread morphogenetic process that filamen- 
tous fungi undertake for dispersal, and also is a key step in the 
colonization of host plants by many fungal pathogens (Roncal et al. 
2002; Kim et al. 2009). Some classical and biochemical studies have 
provided fundamental information on the conidiation in C. carbonum 
(Welz and Leonard 1993, 1994). Nevertheless, little is known about 
the molecular mechanisms controlling the development of conidia in 
this fungus. In the present study, for the first time, we characterized 
a gene BZconl that is critical in controlling the conidiation of C. 
carbonum. BZconl encodes a putative transcription factor (TF) that 
may regulate the conidiation of C. carbonum because we found there 
are two common conserved SANT/Myb-like domains in the N-terminal 
of BZCON1, consistent with the findings from previous studies 
(Grotewold et al. 1991; Wieser and Adams 1995; Aasland et al. 
1996; Boyer et al. 2002). But we still need further functional experi- 
ments to verify this inference. 

Among fungi, in addition to BZconl of C. carbonum, only a Sant/ 
Myb DNA-binding domain-related gene FlbD is known to control 
conidiophore development in A. nidulans (Wieser and Adams 1995), 
but there is no identifiable role for conidiation in the other Sant/Myb- 
like genes (rca-1 of Neurospora crassa and MYT2 of Gibberella zeae) 
(Shen et al. 1998; Lin et al. 2012). MYT2 might have an important 
regulatory role for the genes required for cell proliferation and differ- 
entiation during perithecium development in G. zeae. FlbD is required 
for peridium formation during sexual development but not for asco- 
spore production. In contrast, BZconl mutations have no effect on the 
sexual growth of C. carbonum (Figure 3). 

To date, studies of the molecular mechanism of fungal conidiation 
focused primarily on several model organisms, including the unicel- 
lular yeast, Saccharomyces cerevisiae (Sheu et al. 2000), filamentous 
saprophytic fungi (Aspergillus nidulans and Neurospora crassa) (Park 
and Yu 2012), human fungal pathogens (Candida albicans and Cryp- 
tococcus neoformans) (Saporito-Irwin et al. 1995; Bahn and Sundstrom 
2001; Vallim et al. 2004), and plant fungal pathogens (Ustilago maydis 
and Magnaporthe grisea) (Kruger et al. 2000; Durrenberger et al. 2001; 
Yi et al. 2008). Interestingly, we also found the presence of BZconl 
homologs in Aspergillus nidulans, Neurospora crassa, and Candida 
albicans (Table SI and Figure SI). Thus, these findings suggest that 
C. carbonum should also be a candidate model fungus with great 
potential to provide new insights into the molecular mechanisms in 
the conidiation of fungi. 

Our results showed that BZconl was directly involved in activat- 
ing the development of conidiophores, because the BZconl mutants 
lost the conidia production capacity. In addition, our findings also 



indicated that the presence or absence of 126 continuous nucleotides 
in upstream putative promoter region of the gene BZconl may po- 
tentially mediate BZconl transcription. However, we also noted that 
the strain SI 29 can produce tiny numbers of conidia on PSA (potato 
sugar agar) but not on PDA. Yet, we found that BZconl was not 
expressed in the sexual reproduction of the isolate SI 29. Meanwhile, 
Ml 13-1 (BZconl mutants of Aall3) could not produce any conidia 
on PSA. Thus, it is conceivable that although BZconl is essential for 
the conidiation in C. carbonum, other gene(s) might be involved. 
Future efforts need to be strengthened for discovering more func- 
tional genes modulating the conidiation in C. carbonum. 
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